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ABSTRACT 

In this paper we study in detail the dynamics of flavor transformation for 
neutrinos propagating in the very dense environment of astrophysical compact 
objects as Type II supernova in post collapse phase and proto-neutron stars. The 
analysis is based on the formalism by Strack and Burrows, who introduced the 
generalized Boltzmann equation for Wigner phase space density. In appropri- 
ate limits the formalism reduces to the usual evolution equations for the wave 
functions or for density matrix elements of Liouville equation. We incorporate 
the most important aspects of neutrino propagation physics: the phenomenol- 
ogy of standard oscillations with MSW resonance induced by ordinary matter, 
collective behavior due to self-interaction, which can produce bipolar and syn- 
chronized flavor oscillations, whose relevance is recognized by recent literature, 
the combination of these effects with collisions including scattering, emission and 
absorption of neutrinos. In our numerical simulations a framework with only two 
flavors are adopted together a geometry with spherical and azimuthal symme- 
try. An expansion into Legendre polynomials is adopted. It has the advantage 
on one hand of incorporating the angular dependence, and on the other hand of 
simplifying the treatment of the self-interactions. This method turns out to be 
sufliciently stable for the considered applications and numerically equivalent to 
other developments in term of multipoles series, as discussed in the literature. 

We focus our analysis on the transition from single-angle to multi-angle behav- 
ior of the flavor content in the neutrino outflow. We include both baryon matter 



and electrons and follow the change of the flavor dynamics as their density is 
varied. We distinguish regions close to the neutrino sphere and more distant re- 
gions and in each region we estimate the critical matter density above which the 
neutrino flavor dynamics is suppressed. We argue that the region where the tran- 
sition can occur will move towards the neutrino sphere and the neutrino signal 
will be enhanced in the later stage of the emission epoch. 

PACS : 14.60.Pq , 26.50.+X , 26.60.-c , 95.30.Cq , 97.60.Jd . 



I. INTRODUCTION 



Today it is commonly accepted that the study of neutrinos from supernovae 
with core collapse is an important contribution to development of knowledge 
about elementary particle physics and also to the understanding of the evolu- 
tion mechanism of astrophysical sources [liJei- In fact the supernova in the post 
bounce stage are substantially similar to a black body emitting neutrinos of all 
three flavors, whose propagation in the very dense stellar material is highly sen- 
sitive to parameters such as mass hierarchy or the mixing angles. It may be 
useful in this sense to distinguish two different types of phenomena related to 
the neutrino flavor conversion in terms of distance from the source and discuss 
specifically the effects of conversion to long distance and short distance. The first 
conversion phenomena are caused by nuclear weak interactions of neutrinos with 
the constituents of ordinary matter such as electrons, protons and neutrons via 
charged and neutral currents, even if the typical average value for neutrinos en- 
ergy of about 10 MeV is much less than the mass of the /i and r leptons and only 
the interaction by neutral current takes place. These interactions can be decisive 
in affecting the conversion of the flavor and occur at a a distance from the core 
r > 10^ km M-UJ. The main phenomenon is the celebrated MSW resonance, 
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that was decisive to solve the neutrino puzzle in the solar case [13|,U4. The second 
conversion phenomenon instead occurs at much shorter distances from the core 
of star, r ~ 10 — 150 km, where concentration of produced neutrinos is so high 
to make self-interactions not negligible. Numerical simulations in recent studies 
show that in this situation interesting effects emerge, all characterized by the 
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collective behavior of neutrinos and antineutrinos which are coupled regardless of 



their energy or propagation direction 
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-|22| . These effects are completely differ- 



ent from ordinary oscillations in matter because they involve colliding neutrinos 
with a cross section that depends on the angle between their momenta. These 
particular features introduce a set of coupled non-linear equations and makes pro- 
hibitive an exact solution of the time dependent kinetic equations in the general 



case 
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24| . However reliable approximations have been developed, which have 



lead to a detailed studies of the neutrino flavor dynamics. 



One of the main results 



18H22J of these studies is the prediction that at a 



certain distance from the neutrino sphere a transition from coherent flavor oscil- 
lations to incoherent ones should occur. The coherent oscillations correspond to 
flavour oscillations that are in phase among the different direction along which 
neutrino are propagating, while in the incoherent ones the different directions are 
decoupled and the flavour oscillations are out of phase. Another predicted process 
is a sharp transition from the synchronized oscillations to the so called "bipo- 
lar" ones, that should mainly correspond to the mutual conversion of neutrino- 
antineutrino pairs of one flavour to another one. The two transitions can appear to 
occur simultaneously, and they are the subject of intense theoretical studies. All 
these phenomena can occur at not too large distance from neutrino sphere, since 
otherwise the neutrino density becomes too small, due to the spherical geometry, 
and the momenta of the interacting neutrinos are necessarily focused along the ra- 
dial direction, which reduces the effective neutrino-neutrino interaction strength. 
These processes are strongly affected by the matter physical conditions, besides 
the neutrino and antineutrino flavor compositions. In particular the presence 
of electrons can reduce the effective mixing angle in the matter to a very small 
value. For high enough electron density it is expected that all flavor oscillation 
processes will be blocked by the increase of the electron effective potential on 
the neutrinos due to neutral weak current interaction. Furthermore, closer to the 
neutrino sphere neutrino scattering and absorption can be relevant and modify 
the overall physical picture. Therefore in this case they should be included in the 
theoretical simulations. At larger distance the matter density decrease rapidly 
and both scattering and absorption become negligible, while electrons can still 
play a role. 
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In this paper we analyze the behavior of neutrino flavor content by Boltzmann- 
like equations, following the formulation of Strack and Burrows [25]. We distin- 
guish between the region close to the neutrino sphere, where the role of matter 
can be of decisive relevance, and the region away from the neutrino sphere, where 
the neutrinos flow essentially freely and along straight line trajectories, where 
both neutrino and electron densities are decreasing. We examine the physical 
conditions at different distances from the neutrino sphere and simulate the fla- 
vor evolution under realistic assumptions within the intrinsic uncertainties of the 
problem. 

II. THE FORMALISM 

Different research groups have implemented numerical schemes capable of solv- 
ing the complex equations governing the dynamics of flavor in a general frame- 
work, using different techniques. In particular multi-angle simulations have been 
used, where each trajectory of neutrinos is followed. It has been also explored 
the possibility of reducing the complexity of the calculation using appropriate 
developments of the angular dependence in a series of polynomials, in such a 
way that the final equations are simpler to solve [21]. This is particularly use- 
ful when the initial angular distribution of the neutrinos is highly non-isotropic. 
These studies did not discuss how the phenomenon of flavor modulation due to 
the above mentioned processes is influenced by the effects of neutrino reactions 
in the background matter, i.e. neutrino scattering or emission and absorption. 
In this paper we include in the numerical simulations also the effects of dynami- 
cal decoherence caused by collisions of the neutrinos with the ordinary matter's 
background. The main processes correspond to protons and neutrons that scat- 
ter, absorb and emit particles in addition to electrons that can drastically alter 
the effective flavor mixing angle. The evolution equations that we use in this 



work are the generalized Boltzmann equations of Ref. [25|] that can simulate the 
complex oscillation phenomenology of neutrinos and simultaneously take account 
the neutrino reactions in matter. We solve these equations using the expansion 
in Legendre polynomials and we obtain a suitably truncated set of equations 
that can be numerically handled by the standard software that use an adaptive 
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Runge-Kutta methods. Our main goal is, on one hand, to study neutrino flavor 
evolution near the neutrino sphere, where the presence of collisions and absorp- 
tion, as well as the electron component, could generate decoherence and tends to 
suppress the flavor conversion and oscillations. On the other hand we extend the 
study to larger distance from the neutrino sphere, where neutrinos flow freely but 
the electrons can still suppress the flavor oscillations. 

In the Appendix details are given on the semi-classical Boltzmann-like equa- 
tions of Ref. 25| and their expansion in Legendre polynomials. Here we sketch 
the main elements of the theory and its numerical implementation. The basic 
quantity is the Wigner-Ville distribution in phase space. It is a quasi-probability 
distribution introduced to study quantum corrections to classical statistical me- 
chanics which links the wave function, that appears in Schrodinger's equation to 
a probability distribution in phase space. In second quantization formalism it 
can be defined as 



where ip^ and ij) are creation and annihilation operators. 



Recently |25j a formalism was developed to describe transport processes in 
terms of Wigner functions, which is able to incorporate the quantum phenomenon 
of flavor conversion of neutrinos. The basic equations of the theory can be de- 
rived heuristically merging Boltzmann equation with Heisenberg equation. They 
describe the spatial-time evolution of the matrix elements of the Wigner density. 
Taking the matrix elements in Fock space, considering for simplicity only two 



neutrino flavors, one can write according to Ref. 25 1 



= {ni\p\nj) = I ^'"^ •^'^ I , (2) 
The Boltzmann-like equation can be schematically written 

where the left hand side describes the free streaming of the neutrinos and the 
right hand side includes the processes that involve neutrinos, i.e. the flavor mix- 
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ing matrix, the neutrino-neutrino interaction as well as scattering, absorption and 
emission by matter. The time derivative of the momentum p is caused by the 
gravitational field and can be neglected. Despite this simplification, the complete 
solution of these equations is extremely difficult. Two main approximations can 
be used to obtain from these equations the neutrino fiavor evolution along their 
path away from the star. One can assume that the neutrino emission from the 
star is quasi- stationary. This is a good approximation if the processes that are de- 
scribed by the simulations have a time scale much shorter than the total neutrino 
emission time. In this case one can neglect the time derivative in Eq. and 
the neutrino fiavor evolution is followed along the radial coordinate, assuming 



spherical symmetry. In particular 
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3|-|32| one can assume that scattering, ab- 
sorption and emission are negligible and that neutrinos fiow freely along straight 
line trajectories. Along each one of them the fiavor content evolves due to the 
neutrino-neutrino interaction and the effective mixing angle in matter. 

Another possibility is to approximate locally the time evolution of the dis- 
tribution functions, considering the local matter as homogeneous and in planar 
geometry. This is justified if, within the characteristic time of the considered 
processes, neutrinos can travel for a distance shorter than the length scale of the 
star matter and neutrino density profiles. In this case one considers the time 
derivative and neglect the coordinate derivative of the distribution functions, and 
the equations are solved locally at a given distance from the neutrino sphere. 
For simplicity we will adopt this second scheme, checking that the underlying 
assumptions are indeed well satisfied. However we will implement the method 
by including the increasing focusing of the neutrino fiow at a a given point of 
the radial direction as the distance from the neutrino sphere increases. The solu- 
tion of the equations based on the quasi-stationary assumption within the same 
numerical method of Legendre polynomial expansion is left to a future work. 

If neutrino scattering is negligible, the focusing of the neutrino fiow is dictated 
by simple geometrical considerations [33]. If R is the neutrino sphere radius, 
at a given point at a radial distance r > R the fiow of the emitted neutrino is 
restricted within an angle 6max with respect to the radial direction. According to 
Ref. jssjl, one has 
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COS 9max = \] ^ - (4) 

and each angular direction Op < 6max is related to the angle of emission 6q from 
the neutrino sphere 



m 

where /ip = cos 6p and fiQ = cosOq. The angular distribution of neutrinos must 
be restricted within this angular width at a given radial distance r, see figure 1 
in Ref. 



III. FLAVOR DYNAMICS CLOSE TO THE NEUTRINO SPHERE 

We start our analysis considering the region near the neutrino sphere, where 
the matter density is high enough to produce relevant effects on the evolution of 
neutrino flavor and the neutrino density is the highest one. This region is not 
only relevant by itself but also it can affect the initial conditions for the sub- 
sequent neutrino propagation. Since the size of the neutrino sphere is different 
for different flavors, during the emission process there are regions with an im- 
balance of neutrino contents with respect to the thermodynamic equilibrium. It 
is then of interest to study the flavor evolution under these physical conditions 
and to estimate the characteristic time scale of the possible flavor conversion. 
As physical parameters we consider the total neutrino density p^, = 10^^ cm~^ 
and vary the matter density with typical values in the interval from 10^ to 10^^ 
g/cm^. The electron fraction, if included, is kept to the value 0.4. This is justi- 
fied by microscopic quasi-static calculations of the structure of the after-bounce 
proto-neutron stars. For future discussion Fig. [1] reports a sample of these mat- 
ter density profiles at different entropy of the envelope calculated following the 
scheme of Ref. |3^. The calculation assume neutrino trapping, thermodynamic 
and chemical local equilibrium, and an entropy per particle 5* = 1 in the core 
and 5* = 1,2,3,4 in the envelope. The lowest entropy profiles should correspond 
to the latest stage of the neutrino emission period. The results are in line with 



dynamical simulations of the Livermore group, see Fig. 1 of Ref. 351] ■ Despite 
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the quasi-static calculations do not extend to the more far tails of the profiles 
and the two sets of calculations use two different types of scheme, they show a 
very similar trend and their overall evolution suggests the range of density values 
that the profile should cover during the neutrino emission epoch. 

For simplicity we consider a two-flavor framework, with parameters corre- 
sponding approximately to the electron and atmospheric neutrinos, see Table [H 



conventionally indicated as electron and muon neutrinos 



For absorption and 



emission processes one needs to specify the neutrino temperatures. We use the in- 



dicative values reported in Table HI taken from Ref. 26|. First we consider equal 
number of neutrinos and anti-neutrinos, but with an excess of electron flavor. 
Neutrinos are considered mono-energetic, and the equations are correspondingly 
adapted as detailed in the Appendix. Following the scheme described in the pre- 
vious section, in Fig. |2] is reported the time evolution of the fractional content 
in electron flavor of the neutrinos, with energy E = 10 MeV. This content P,y^ 
is normalized to the total number of neutrinos and anti-neutrinos. The angular 
distribution at the neutrino sphere is taken "half isotropic" , i.e. uniform angular 
distribution in the forward direction, and no backward emission, to simulate the 
flow from the neutrino sphere. We take the radius of neutrino sphere equal to 
10 Km and analyze the region around a radial distance of 15 Km. We start the 
simulation from an angular distribution "focused" according to Eq. (jl]). The aim 
of the calculation is to check the stability of the angular distribution with time, 
in order to estimate the time necessary for the appearance of the transition to 
the multi-angle regime. If the matter electron component is neglected as well as 
the emission and absorption processes, the temporal evolution of P^^ is depicted 
in Fig. [2^. The time unit is taken to be the oscillation time period in vacuum 
Ty, as specified in the Appendix. In a time as short as 10~^ s there is a sharp 
drop. It corresponds indeed to the transition from the "single-angle" evolution 
to the "multi-angle" one. This transition is apparent in the three-dimensional 
plot of Fig. [3], where the angular distribution of the electron flavor content as a 
function of time is reported. Initially the angular distribution stays stable with 
time. At the transition time, apparent in the upper panel of Fig. |3l the angular 
distribution changes abruptly and becomes irregular and time dependent. It is 
to be noticed, however, that the angular distribution after the transition is not 



completely random, but displays some structure. Furthermore it has to be kept 
in mind that the angular distribution is a probability density and it can acquire 
arbitrary values. Notice also the short time necessary for the transition, which 
justifies the assumption of uniform matter that we are using. 

If we introduce the electron component at increasing density, the signal dis- 
appears at some critical value. This is expected since the effective mixing angle 
tends to vanishing small values as the electron density increases. Correspond- 
ingly, the angular distribution remains unchanged with time, as depicted in the 
lower part of Fig. |31 The actual matter and electron density depends on the stage 
of the neutrino emission epoch. At a fixed distance in the initial stage the matter 
density is likely to be larger than the critical one, while at a later stage it can be 
smaller and the neutrino signal should then appears. 

In the second column of Fig. [2] the same analysis is reported for the inverted 
hierarchy for the neutrino mass spectrum. A similar behavior is observed, but 
the transition occurs at an early time. 

The same type of calculations of Fig. |2] are reported in Fig. H] for a neutrino 
density py = 10^^ cm~^. The observed trend is quite similar. However the tran- 
sition is obtained at a larger matter density and the time scale for the transition 
is longer. The time scale for the occurrence of the transition is still acceptable for 
the validity of the method of analysis. One notices the scaling behavior : as the 
density decreases by two orders of magnitude, both the transition density and the 
time for the transition increase by one order of magnitude. The corresponding 
angular distributions are reported in Fig. [51 Also in this case the appearance of 
the transition is quite evident. 

In all these cases we found that both absorption/emission and scattering pro- 
cesses have essentially no relevance. To put in evidence their possible role, we 
suppressed the electron component and we repeated the calculations at the same 
neutrino density but at higher matter density and at two different neutrino ener- 
gies, assuming direct hierarchy. The results can be seen in Fig. |6l They are able 
to suppress any flavor dynamics, but only if the matter density is much higher. 
As a consequence the inclusion of the electron component, with the correspond- 
ing density, would wash out any signal, as we have explicitly checked. Therefore 
the absorption/emission processes are expected to have very limited effects in all 
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cases. The scattering processes have a much higher time scale and play no role in 
any case. To illustrate further this point we consider the evolution of the flavor 
content for an isotropic angular distribution, which would be more appropriate 
for the neutrinos inside the neutrino spheres of both flavors. Accordingly, the 
matter density is taken at 10^^ g/cm^ and the neutrino density at 10^^ cm~^ . 
We assume the same initial flavor contents as in the previous calculations, so that 
the simulations should give the characteristic time scale for the relaxation toward 
equilibrium. For aim of comparison, we take the same temperature values for the 
two flavor neutrinos as in the previous analysis. As reference, vacuum oscilla- 
tions are reported in Fig. [7^. If only the neutrino self-interaction is included the 
usual bipolar oscillations are observed, see Fig. [TJo. If only absorption/emission 
are included, the time scale of the flavor dynamics increases by three order of 
magnitude, see Fig. In this case one gets a relaxation toward equilibrium 
through damped oscillations. Despite the time scale is drastically enhanced, it 
is still much shorter than the ones of typical macroscopic processes inside the 
neutrino sphere, as the trapping time or the possible density oscillations. If in 
addition also the self-interaction is introduced, apparently no change is observed, 
i.e. the absorption/emission processes dominate the flavor evolution, see Fig. [T]!. 
Finally the inclusion of the electron component damps all the oscillations in the 
evolution toward equilibrium with or without self-interaction. Figs. Uh, f. This 
analysis, even if only illustrative, indicates that the time scale toward flavor equi- 
librium is much longer than typical bipolar oscillation period but shorter than 
the time scale of the macroscopic processes inside the neutrino sphere. 

We consider now the effect of the neutrino-antineutrino asymmetry. Accord- 
ing to Ref. {321 if the flavor content is substantially different for neutrino and 
antineutrino, the flavor dynamics and the transition from single-angle to multi- 
angle regime are suppressed. Assuming initially zero flavor content of /i neutrino, 
i.e. Pui^ = Pu^ = at t = 0, the asymmetry can be specified by the value of the 
parameter e defined as 

P - P- 

equal to the relative fraction of different neutrino and antineutrino electron flavor 
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content. 

We have investigated this effect under the physical conditions of Fig. [7^. The 
variation of the flavor time evolution at different asymmetries is reported in Fig. [HI 
The suppression of decoherence seems to start around e = 10'^. We therefore 
confirm the phenomenon of decoherence suppression by asymmetry. The precise 
value of the critical value of e depends of course on the effective mixing angle and 
on the neutrino density. In general a value of the asymmetry as small as 10"^ or 
so is enough to fully suppress the flavor coherence. 

IV. MOVING AWAY FROM THE NEUTRINO SPHERE 

We will consider in this section the flavor dynamics at larger distance from 
the neutrino sphere. With respect to the previous analysis the main change that 
occurs is the degree of the focusing of the interacting neutrino beam. We do not 
report all the results, but we concentrate on the case of neutrino density 10^^ 
cm~^. We consider the distance of 50 km, away from the neutrino sphere. One 
can see from Fig. [9] that the time evolution at different matter density of P^^ is 
more irregular. Despite that, it turns out that the transition is still present, as it 
can be seen from the angular distributions of Fig. [TUl corresponding to Fig. [Hfe, 
i.e. at a matter density 10^ g/cm^. The transition is much smoother and it is 
not clear if a real multi-angle regime is reached. For the inverted hierarchy the 
evolution is more regular. To be noticed that at a density slightly smaller than 
the critical one, where the signal disappears, the oscillations are quite regular. 

It can be interesting to look at the time evolution of the angular distribution 
in the case where no matter is included, i.e. for a gas of free neutrinos with self- 
interaction only. The time evolution of the angular distributions corresponding 
to 30 and 40 km. are reported in Fig. [TTl In this case the only effect is coming 
from the focalization due to the distance. One can observe that there is still a 
transition of regime, but it is much smoother in time than in the previous cases. 
The angular distribution changes smoothly from a uniform one to a structured 
one. The latter, however, does not appear so random as in the previous cases 
but it displays only some modulations. The mixing angle has the vacuum value, 
much larger than the one in the matter with electrons, and this can be the reason 
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of this different behavior. 



V. DISCUSSION AND CONCLUSION 

For a given flavor composition and symmetric neutrino-antineutrino flavor con- 
tent, we have estimated the critical matter density below which the transition 
from single-angle to multi-angle behavior can be expected and the time scale 
of the transition for a wide set of density values. The method is based on the 
expansion of the Boltzmann equations in Legendre polynomials. The calcula- 
tions were performed both near the neutrino sphere and away from the neutrino 
sphere. The focusing of the neutrino flow at different distances modifles both the 
critical density and the time scale. Once a density proflle is given, one could get 
an estimate of the region where the transition should flrst occur. Since the den- 
sity proflle changes with neutrino emission time, this region will move with time. 
In general the matter density proflle is uniformly decreasing along the emission 
time, therefore the region for the transition will move toward the neutrino sphere. 
Furthermore since the neutrino flavor modulations are more apparent near the 
neutrino sphere, the neutrino signal should be stronger at the later stage of the 
neutrino emission epoch. The transition appears to be smoother as the distance 
increases, due to the focalization of the neutrino beam. If the neutrino flow is 
asymmetric, the flavor evolution changes drastically as the asymmetry reaches a 
critical value, that is estimated to be not larger than 0.1. This is in agreement 
with the flndings in Ref. 32|. This type of analysis is usually performed by in- 
troducing the so called "polarization vector" , whose three-dimensional evolution 
in time (or distance) gives an overall view of the flavor dynamical evolution. We 
preferred to use directly the full angular distribution of the flavor content since 
it gives a more detailed description of the flavor structure. Of course the two 
methods are equivalent and correspond to two different representations. 

The analysis was done within the approximation of local uniform matter for 



the Boltzmann equation, where the time is the only evolution parameter 
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20|. If a transition is present, it occurs in a very short time, which justifles 
the approximation. However another approximation method is the stationary 
assumption, where the only evolution parameter is the distance from the neutrino 
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sphere 
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25 



32 



33| . This is justified since the time evolution of the density 



profile is much longer than the neutrino transit time through the envelope. This 
alternative analysis within the polynomial expansion, and its comparison to the 
local approximation, is left to a future work. 
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Appendix: Transport with all main neutrino interactions, symmetry con- 
ditions and equations on the set of Legendre polynomials 

Recently it was developed a formalism to describe transport processes in terms 
of Wigner functions, which is able to incorporate the quantum phenomenon of 
flavor conversion of neutrinos [37|. In terms of the matrix of equation (|2]) the 
complete evolution equation is 

In this equation diagonal terms ere real valued quantities that have the mean- 
ing of phase space densities, while off-diagonal terms are complex quantities with 
the meaning of macroscopic overlap functions. is the mixing Hamiltonian 
and C is the collision matrix. In general the Hamiltonian includes three differ- 
ent terms corresponding to vacuum mixing, that is energy dependent, to matter 
interaction, that depends on the density proflle, and the neutrino-neutrino self- 
interaction. The last contribution is more difficult to treat because it depends 
on the angle among neutrinos impulses. Introducing real and imaginary part of 
off diagonal macroscopic overlap as fr and /j, one can easily find the generalized 
Boltzmann equations for two neutrino fiavor interacting with a background and 
with neutrino-neutrino interactions included, and similarly for their antiparticles. 

For application it is useful rewrite transport equations in terms of specific in- 
tensities of neutrino radiation field using the relation connecting invariant density 
fy to corresponding specific intensity I^, according to ly = E^f^/ {2Txfif c^. 
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The main terms to be discussed in detail for our purposes are certainly collision 
term and neutrino- neutrino interaction term By. 

Generic collision term Cy is usually written in terms of specific intensity. For 



convenience we report here the explicit expression taken from Ref. 
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where = {2-Khf CyC/ . 

The quantity is a phase function for scattering integrated over the solid 
angle and gives the probability that a particle enters into the beam or go out. 
For a scattering process i it is well approximated by 



(fi, Vt') = ■ Vt') = I + 6iVt ■ Vt' = I + di^iy, (A.3) 

where 5i is a constant specific to each scattering process and /ijy is the angle 
between the incident and outgoing neutrinos. The quantity is the sum of all 
absorption processes, rii is the number density of matter species i and erf denotes 
the absorption cross sections. Similarly for scattering processes one introduces the 
corresponding quantities and af. The quantity J-'^'' is the equilibrium Fermi- 
Dirac occupation probability for neutrino species and is the corresponding 
black body specific intensity 



Generic neutrino self-interaction term By has an integral form and using the 
distribution fy it can be written as 

By (p, v,t) = f3 j d\ (1 - COS e^^) fy (q, r, t) , (A.4) 

where 6^*^ is angle between a test neutrino with momentum p and a fixed back- 
ground neutrino with momentum q. In this coupling coefficient integration goes 
over all momenta in the ensemble of neutrino background and the strength is 
given 

Purely kinematic aspects about neutrinos transport are in the partial deriva- 
tives at first member of the discussed equations while at second member the 
dynamical terms governing fiavor conversion and collisions are present. In six- 
dimensional phase space the single particle Boltzmann equation can be generally 
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expressed in compact form as the total derivative for an arbitrary / equal to 
collision integral. 

We consider the simplest case in which neutrinos are distributed uniformly 

in space and can be scattered by the nucleons when simultaneously interact 
with themselves and electrons. We also restrict neutrino energy spectrum to 
be monochromatic. 

It should be noted that restricting the energy to a single value implies a choice 
for the structure of the emission term. Indeed, this term depends on the neutrino 
energy and the temperature of a specific flavor, but is not proportional to the 
distribution function. Assuming the / proportional to a Dirac delta function in 
energy to select a fixed value in the spectrum a choice can be made for the emission 
term and two different types of approximation are possible. The first implies an 
integrated emission term not energy dependent, while the latter implies energy 
dependence via cross section. In this paper we fixed temperature to typical values 
and use the second choice. 

With all discussed approximation the time derivative of the distributions / is 
only term of the left-hand sides of equations not vanishing in a fixed cartesian 
system of coordinates. 

The mathematical structure of transport equations and the angular depen- 
dence of the densities provides an useful expansion on the set of orthogonal Leg- 
endre polynomials. In fact different terms depend on the unknown functions with 
a kind of linear or quadratic dependence. Specifically the self- interaction neutrino 
term has integral form and multiplies other terms containing the functions /, 
which results in quadratic expressions. A generic fi, for a fixed energy value has 
components defined as the integrals 

/» {t, E) = d//pP« (/xp)/. (t, E, /xp) , (A.5) 
and it is possible to reconstruct the original densities using the expansion 

U {t, E, = E + ^^'^ M /i'^ ^) ■ 

The typical time scale characterizing neutrino self-interactions is the period 
associated with bipolar oscillations, tb — 27r/^/TJ]2, where /i — ^rii, with rii, 
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the numerical density of all neutrinos species. It is also useful to introduce the 
period of free oscillation in vacuum Ty — L/cin terms of oscillation length L — 
AirhE/ Am^c^ depending on energy and neutrino mass square difference. For 
convenience we choose Ty as unit of time, and define a specific interaction rate 
due to collisions in terms of times r/ = 1/cNApYia^ for scattering and rf = 
{1 — J^l'^) / cN ApYiCT^ for stimulated absorption-emission. The quantity Yi are 
fractions of nucleons, o"? and cr" are cross sections. Index i indicates nucleon and 
can be equal to n for neutrons, p for protons and t for the sum of processes, 
with T^~^ — Tp~^ + T^~^. Note that interaction rates are not flavor independent 
because the cross sections depend on the mass of the lepton specifically involved 
in the process. 

With above assumptions, if all distribution and overlap functions are normal- 
ized to the neutrinos density, the equations can be rewritten in a form useful for 
their numerical solution. For neutrinos and antineutrinos one gets a set of eight 
coupled nonlinear equations. For example, the first of these has the explicit form 

/iP = 27ri^(')sin2^+ 
+2n ((/(°) - F) /«- (/f - F) f^' + 

(A.7) 

- (f ) /z^P + (^) {b* Si,o - A']) + 

+ (f ) {f^%,o + + (^) {fVs,o + 'ff^%,.). 

The angle 6 is the vacuum mixing angle between two neutrino species and the 
density of neutrinos appears due to the normalization. The quantity B* {T^J is 
the integrated emission term depending on temperature according to the above 
discussion. The choice of time unit is convenient since the ratio tb/tv is depen- 
dent only on the ratio ui/ p. 
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In four of the eight equations a parameter A appears, that represents the 
potential induced by ordinary matter. It is proportional to concentration of 
charge lepton according to A = 2\/2LGFne/ {nhc). Note that above equations 
have some terms with different algebraic signs from the original reference. This 
depends on different convention used for the relative sign between the real and 
imaginary parts of complex functions associated with neutrinos and antineutrinos. 
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Parameter 


Value 




7.59 • 10^^ eV^ 


sin (29) 


0.93 




2.76 MeV 




6.26 MeV 




4.01 MeV 


T- 


6.26 MeV 




-0.2 


K 


-0.1 



TABLE I: Neutrino physical parameters used in the simulations. For the meaning 
the parameters 5 see Eq. (lA.Sp . 
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5 10 15 20 25 30 

r{km) 

FIG. 1: Baryon density profiles of a proton-neutron star of mass M/Mq = 1.4 for 
different entropy. The core entropy is fixed at S = 1 while the envelope entropy is 
taken, from the lower to the higher profiles, asS = l,2,3,4,5. 
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FIG. 2: Time evolution of Ve fraction at radial distance of 15 km and different densities 
matter values. Energy is fixed to 10 MeV and neutrino density to 10'^^ cm^'^. Time 
is in unit of vacuum oscillation period Ty- Left: normal hierarchy. Right: inverted 
hierarchy. 
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0.010 0.0 




/v. 



0.010 0.0 



FIG. 3: Angular distribution of Vf. content in function of time with same parameter of 
Fig. [2] in normal hierarchy. Quantity f^^ is normalized to total neutrino density and 
/ip = cos(6'p), with 6'p the angle with respect to the radial direction Upper plot: density 
matter fixed to 10'' gr/cm'^. Lower plot: density matter fixed to 10^ gr/cm^. 



22 




0.000 0.002 0.004 0.006 0.008 0,010 0.000 



0.004 0.006 OOOK 0.010 



0.00 0.02 



p = 1 . X 1 0"'' gr/cm.^ 




0.06 0.08 



1,0 -(d) 



0.10 0,00 



p - l.x 10-^ gr/cm^ : 




0,04 0,06 0,08 0.10 



10 ■ (e) 



p= l.xlO'-gr/cm' : ,„:(f) 




p = 1.x lO' gr/cm' : 




10 0.00 



0.04 0.06 0.08 0.10 



10 Kg) 



p= l.x lO' gr/cm' 




-(h) 


p= l.x lO' gr/cm' . 







0.000 0.002 0.004 0,006 



0.010 0.000 0.002 0.004 0.006 0.008 0.010 



1.0 ■ (i) 



0.000 0.002 



0.004 0.006 
t/ry (Ty = 1.09x10"' s) 



0.008 0.010 0.000 



p= 1 . X 1 gr/cm 



0.004 0.006 0.008 0.010 

t/Ty lTy= 1.09x10"^ S) 



FIG. 4: The same as in Fig. [2] but for neutrino density fixed to 10^^ cm ^. 
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0.010 0.0 




0.010 0.0 



FIG. 5: Angular distribution of fe content as a function of time with the same parame- 
ters of Fig. m in normal hierarchy. Quantity f^^ is normalized to total neutrino density. 
Upper plot: density matter fixed to 10^ gr/cm'^. Lower plot: density matter fixed to 
10® gr/cm'^. 
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(b) E =15 MeV, =0, p = 1.x lO' gr/cm^ - 




0.0000 0.0001 0.0002 0.0003 0.0004 0.0005 0.0000 0.0001 0.0002 0.000.1 0.0004 0005 
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FIG. 6: The same as in Fig. [2] but at higher matter density and for two different energies 
in direct hierarchy. Neutron fraction is fixed to 1 to put in evidence possible role of 
absorption, emission and scattering. Left: Energy is fixed to 10 MeV. Right: Energy 
is fixed to 15 MeV. 
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FIG. 7: Time evolution of fraction for initial isotropic angular distribution and 
different physical conditions. Energy is fixed to 5 MeV, neutrino density to 10^^ cm~^, 
density matter to 10^^ gr/cm^ and hierarchy is direct. Time is in Ty unit. Panel 

(a): vacuum oscillation, (b) self-interaction only, (c) only absorbtion/cmission, (d) 
self-interaction and absorbtion/emission, (e,f) including also electrons with or without 
self-interaction. 
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FIG. 8: Neutrino flavor dynamics at different values of the asymmetry parameter e, see 
the text for details. 
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0.90 

FIG. 10: Angular distribution of Vf, content as a function of time with the same pa- 
rameters of Fig. [9^. 
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FIG. 11: Angular distribution of Vf, content as a function of time for a free neutrino 
gas with only self-interaction at the distance of 30 km (upper plot) and 40 km (lower 
plot). Neutrino density is fixed to 10^^ cm~^ and hierarchy is normal. 
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